Continuity of the Explosive Percolation Transition 
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The explosive percolation problem on the complete graph is investigated via extensive numerical 
simulations. We obtain the cluster-size distribution at the moment when the cluster size hetero- 
geneity becomes maximum. The distribution is found to be well described by the power-law form 
with the decay exponent r = 2.06(2), followed by a hump. We then use the finite-size scaling 
method to make all the distributions at various system sizes up to A*' = 2"^^ collapse perfectly onto a 
scaling curve characterized solely by the single exponent r. We also observe that the instant of that 
collapse converges to a well-defined percolation threshold from below as N ^ <x. Based on these 
observations, we show that the explosive percolation transition in the model should be continuous, 
contrary to the widely-spread belief of its discontinuity. 

PACS numbers: 64.60.ah, 64.60.aq, 36.40.Ei 



The term explosive percolation was proposed in Ref. [1] 
to describe a sudden appearance of a macroscopic clus- 
ter in a network growth model with the so-called product 
rule considered on the complete graph. This growth rule, 
named as the Achlioptas process (AP), is then studied 
on the two-dimensional lattice 0, Q and on the scale- 
free networks [4-6] as well, yielding similar results. That 
suddenness has been widely believed to indicate a dis- 
continuity at the percolation transition in the thermody- 
namic limit 0, Q : and the similar explosiveness has been 
observed with the other growth rules proposed later 
[isj . These observations of the explosiveness have drawn 
much interest due to the striking difference from the well- 
known continuous transition in the standard percolation 
models ^14] . However, in our point of view, the explosive- 
ness has not been carefully investigated as yet enough 
to draw a decisive conclusion on the discontinuity, and 
possibly just represents an extremely steep but still con- 
tinuous transition. 

Friedman and Landsberg Q have suggested the argu- 
ment of the powder keg as a circumstantial description 
to explain the apparent discontinuity of the explosive 
percolation transition. Meanwhile, da Costa et al. (isj 
have reported that the explosive percolation is actually 
continuous for a modified version of the AP by analyti- 
cally deriving the critical scaling relations based on nu- 
merical observations of power-law critical distribution of 
cluster size [l6| . In this Letter, we try to unmask the 
(dis)continuity in a systematic and direct way by per- 
forming a careful finite-size-scaling analysis at newly in- 
troduced psewrfo-transition points for finite systems and 
show that the explosive percolation transition on the 
complete graph is indeed continuous in the thermody- 
namic limit. 

The model we study is_the AP with the product rule 
on the complete graph Start with N nodes with 
all links unoccupied. At each step, choose two possible 



unoccupied links randomly between nodes. Then, select 
the link merging two clusters with a smaller product of 
the two cluster sizes. Here, a cluster is defined as a set 
of nodes connected each other via occupied links. This 
procedure is repeated until all nodes are connected as a 
whole. The number of occupied links L increases one by 
one at each step and the occupied link density (or time), 
t = L/N, serves as the control parameter for the model. 
The interested observable is the largest cluster size G{t) 
which becomes macroscopic (linear in N) at sufficiently 
large t. The order parameter is defined as the relative size 
of the largest cluster, g{t) = G{t)/N, which remains at 
zero below the threshold tc and becomes finite for t > tc 
in the N = oo limit. 

The main question is whether the gap, g{t)\^_^^+, van- 
ishes (continuous transition) or approaches a non-zero 
constant (finite jump). It may be natural to use the in- 
formation above the transition point (< > tc) in order 
to prove the (non-)existence of the gap or estimate the 
gap size. Thus, most of previous studies have focussed 
on this information [l|-[l3l [is} , but could not provide a 
definitive answer due to the extremely slow convergence 
of the order parameter in system size. In this work, we 
took the opposite approach. Using the accurate infor- 
mation below tc, it is still possible to derive the upper 
bound for the gap, which turns out to vanish as iV ^ oo. 
This guarantees the vanishing gap at the transition. 

Our strategy is as follows: (i) Set up lower and up- 
per pseudo-transition points, ti{N) and t„(7V), for finite 
size N below and above the true asymptotic percolation 
transition point respectively. We expect that both 
pseudo-transition points converge to as TV — oo. (ii) 
Find the upper bound for the size increase of a largest 
cluster AG between ti{N) and tu{N). (iii) Show that 
this upper bound is sublinear in N, which implies no 
macroscopic jump of the largest cluster size through the 
percolation transition. This completes the proof of the 
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FIG. 1: (Color online) Convergence of ti (the lower branch) 
and tu (the upper one) averaged for 100 ~ 5000 different re- 
alizations for each system size (A'' = 2^'^ ,2^^ , ..,2^''). Errors 
are smaller than symbol sizes, if not shown explicitly. Lines 
are just guides to the eyes. As A'' increases, tu (ti) decreases 
(increases), approaching to the well-defined value in the ther- 
modynamic limit displayed as the horizontal line. For clarity, 
the upper branch is vertically enlarged in the inset. 



continuity at the explosive percolation transition. All 
procedures are done via extensive numerical simulations, 
typically up to = 2^"^ ss 1.37 x 10^-^ and the average is 
done over 100 ~ 5000 runs. 

The most crucial step is to define the two pseudo- 
transition points at the microscopic step level. First, we 
introduce the lower pseudo-transition point ti (N) as the 
instant when the cluster size heterogeneity (the number 
of distinct cluster sizes) becomes maximum. For small t, 
the cluster size heterogeneity increases with t due to the 
randomness of clustering processes. However, the emer- 
gence of a macroscopic percolating cluster which contin- 
uously absorbs small clusters causes the heterogeneity to 
decrease and eventually the whole system becomes one 
cluster. Due to the mechanism of suppressing the emer- 
gence of large clusters, one may argue that the hetero- 
geneity increases slowly but steadily up to just before 
the explosion when many different size clusters merge 
into one big macroscopic cluster. Thus it is reasonable 
to consider the maximum heterogeneity as a preceding 
symptom of the percolating onset for finite systems. In 
Fig.[l] the average values of ti{N) are plotted against N 
(lower branch), which converge to the asymptotic value 
of tc = 0.8884490(5) from below, as expected. 

Second, we expect that the growth rate of the largest 
cluster also becomes maximum at the percolation transi- 
tion. Microscopically, the upper pseudo-transition point 
tu{N) is defined as one step after the moment when the 
second-largest cluster size becomes maximum. Thus G(t) 
can experience a largest increase exactly when t exceeds 
tu, since the second-largest cluster merges into the largest 
cluster. A typical growing process is displayed in Fig. [21 
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FIG. 2: (Color online) Evolutions of the largest cluster size 
(the upper curve) and the size of the second largest cluster 
(the lower curve) versus the growth step t. As t crosses tu from 
below, the size of the largest cluster exhibits a sudden biggest 
increase since the maximum second-largest cluster merges into 
it at that moment. The system size is A'' = 2^'^. 



Note that the second-largest cluster never recovers its 
size after merging into the largest cluster. So there will 
be no explosive increase of the largest cluster size for 
t > tu{N). Dominance of one percolating cluster is the 
characteristic of the percolating phase. So it is reasonable 
to expect that tu{N) is just above tc, which is consistent 
with numerical results (see the inset in Fig. [T]). The av- 
erage values of tu{N) converge to the same asymptotic 
value of tc from above as A^ oo. 

The sample-to-sample fluctuations decrease with 
^-0.5 j-j-^Q^ shown here), which implies that both i„ and 
ti are self-averaging [l7|, so not only the critical point 
but also any sample-averaged quantity are well defined 
in the asymptotic limit [16]. We also find numerically 



tu{N) - tiiN) N- 



(1) 



with 5 = 0.39(3) [181. 

Now we consider Ag = g[tu) — g{ti) — AG/N, which 
is the growth of the largest cluster density between two 
pseudo-transition points through the asymptotic transi- 
tion point tc- In the thermodynamic limit, it will be the 
jump size (if any) of the order parameter at the percola- 
tion transition. As g{ti) must vanish as A'^ ^ oo, we only 
need the information of g{tu) in principle to calculate 
A^Iat^oo- Figure [3] shows the cluster-size distribution 
n(s; tu) of cluster size s (normalized by the total number 
of clusters C) at tu{N) for various different sizes A^. The 
distribution fits extremely well with a power-law form, 
n{s;tu) ^ s~'^ with the decay exponent t — 2.06(2) in a 
huge range, which is then accompanied by a little dip in 
the end. In Fig. O the largest cluster size G{tu) depend- 
ing on the dip structure near the upper cutoff shows a 
slight trend of the sublinearity in N (moving left in the 
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FIG. 3: (Color online) Log-binned cluster-size distributions 
at tu for each A'", where the horizontal axis is the cluster size 
s divided by N for convenience. The solid line above is a 
guiding line of which the decay exponent is 2.06. 



axis of s/N as N increases), which may be one symptom 
for the continuous transition. However, as discussed be- 
fore, it can not be conclusive even with huge system sizes 
studied here. 

If one assumes a conventional natural cutoff of the 
power-law type distribution function, the upper cutoff 
which should be proportional to the largest cluster size 
G{tu) will scale as N^/(-'-^') ~ N°-^'^ with t ~ 2.06. Sub- 
linearity is estimated only by 6%, which may call for a 
huge system size like N 10^'' (beyond the present com- 
puting capability) to reach a reasonable scaling regime 
{g{tu) ^ 0.1) and get any sensible extrapolation to the 
thermodynamic limit. Most of previous studies [l- 13. 15| 
basically depend on the data in this supercritical regime 
(t > tc). Nevertheless, the scaling plot with this natu- 
ral cutoff shows a reasonable collapse including the dip 
structure at the end, but involving big statistical errors 
(not shown here). 

In efforts to find conclusive evidences, we scrutinize 
the cluster-size distribution n(s; ti) at ti(N), which shows 
again the power-law decay with the same decay exponent 
T = 2.06(2) followed by a hump near the upper cutoff (see 
Fig. |4]). It has a much shorter (but still quite broad) 
power-law regime, but exhibits much better statistics 
even in the hump region (see Fig [5]) . This power-law 
scaling with a hump-like structure at the end has been 
also reported in previous studies [1, H, H, [13, El] • 

In contrast to the cluster distribution at i„(7V), n(s; ti) 
shows a fast exponential decay near the cutoff s / . This 
sharp cutoff originates from the nature of the growth 
(product) rule which discourages the merging of bigger 
clusters before explosion. One may estimate Sf{N) at ti 
as follows. It is appropriate to estimate the upper cutoff 
s/ by assuming the 0(1) number of clusters left beyond 
the cutoff, i.e. J2s>sf '^('*) ~ where C is the total 



FIG. 4: (Color online) Log-binned cluster-size distributions 
at ti for each A^. The decay exponent of the guiding solid line 
below is 2.06. 



number of clusters in the system. Note that C scales lin- 
early with N |19[. Since n(s) decays exponentially fast 
(or faster) near s/, X]s>sf "'('^) ~ (^'S/)'^(s/) with a fi- 
nite characteristic scale Asf of the fast decaying part. 
Consequently, one can find n{sf) ^ sj'^ ^ 1/C ^ '^/N, 
which leads to 



ss{N)^N^'\ 



(2) 



at t = ti {N). Note that this cutoff scales differently from 
the natural cutoff. The huge difference in the largest clus- 
ter size G{t) just below and above tc leads to its abrupt 
and explosive increase through the percolation transition, 
which is the main difference between the ordinary and 
explosive percolation. However, the magnitude of the 
explosion may be still sublinear in N as discussed before. 

The validity of Eq. ([2]) can be checked numerically 
through the finite-size-scaling (FSS) analysis for the en- 
tire distribution function n{s; ti). Assuming a single char- 
acteristic cluster size (proportional to Sf), one can write 
a FSS form for the distribution n(s]ti) as 



n{s;ti) = s--f{s/sf) = s'^ f{sN-'/-), 



(3) 



where f{x) is the scaling function that becomes an 0(1) 
constant for a; ^ 1 and decays exponentially fast (or 
faster) for x > 1. In Fig. [5l we plot s'^n{s;ti) averaged 
over 100 different runs versus the scaling variable sN~^/'^ 
for aU 21 different system sizes N = 2^'^, 2^^, • • • , 2^^, 2^'^, 
using T — 2.06. As is clearly seen, the collapse of the 
data points is perfect including both the power-law decay 
part and also the hump structure near the end. This re- 
markable collapse validates the single- variable FSS form 
of Eq. ([3]) without any doubt. Therefore, we now have 
the most precise and full information on the cluster size 
distribution just below the transition for large N . Com- 
parison with the scaling collapse plot using r = 2 (see 
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FIG. 5: (Color online) Finite-size scaling of the cluster size 
distribution n(s) at ti{N): All data points plotted in the form 
of s^n(s) versus sN^^^^ at r = 2.06 collapse into a single 
curve. Inset: The same plot but with r = 2. This clearly 
shows that r ~ 2.06 (> 2) gives us the much better quality 
of the FSS collapse. 



the inset of Fig. [5]) leads to the definite conclusion that 
T must be larger than 2. 

Now we are ready to derive the upper bound of 
^5 = gi^u) — g{U)- From t = ti to t = tu, we need 
AL = N{tu — ti) steps {AL links added). One may imag- 
ine the ideal process to maximize the growth of the largest 
cluster G{t), starting from the well-known cluster distri- 
bution rt(s; ti) at t;, by adding AL links one by one. This 
ideal process can be easily implemented by simply link- 
ing and merging the largest cluster with the next largest 
cluster at each step and repeating it till all AL links are 
exhausted. Then, all clusters of size s > sg will merge 
into one cluster, which becomes the largest cluster af- 
ter AL steps. The threshold value sg is determined by 
balancing the total number of merged clusters with the 
total number of links added; C X]s>s5 "^(^i ~ with 
C the total number of clusters at t = t; . 

During this ideal process, the largest cluster G(t) grows 
by the amount of C^^^^^ sn{s;ti). One can easily esti- 
mate ss ^ (C/AL)i/(i-^) ~ N^/iT-i)^ using the single- 
variable FSS form of n{s; U) of Eq. Q with Eq. Fi- 
nally we get the strict upper bound for Ag as 



A5 



< 



(4) 



This shows that the order parameter jump Ag at the 
percolation transition should vanish as the N ^ oo limit, 
if (5 > and t > 2, which are undoubtedly confirmed in 
our numerical simulations. Therefore, we conclude that 
the explosive percolation transition is indeed continuous. 

In summary, we showed that the explosive percolation 
transition on the complete graph is continuous by ex- 
ploiting the high-precision cluster-size information at the 
moment of the maximum cluster heterogeneity, i/(iV), 



approaching the asymptotic transition point tc from be- 
low. The cluster-size distribution displays the power-law 
scaling with the decay exponent r — 2.06(2), followed 
by a hump with a sharp cutoff s/ N^/'^. It is ex- 
plicitly shown that the existence of the single-variable 
finite-size scaling at i/ {N) solely guarantees the continu- 
ity of the transition if r > 2. Therefore, the scaling and 
the discontinuity can not be compatible near the explo- 
sive percolation transition as in usual critical phenomena. 
The explosiveness originates from the huge difference in 
the largest-cluster-size scaling in TV below and above the 
transition. However, it is not enough to invoke a discon- 
tinuity at the transition. 

Our approach can be applied to many other models 
including various different types of explosive percolation 
models to clarify the (dis)continuity. Applications to 
other explosive percolation problems and also the low- 
dimensional cases are currently under investigation. 

Note added. A few days before our submission, Grass- 
berger et al. posted a preprint [20] drawing a similar con- 
clusion, but following a completely different approach. 
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